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I demonstrate that instability in a system of homogeneous scalar fields leads to the growth of 
super-Hubble metric perturbations. This generalizes the result that parametric resonance can lead 
to the growth of cosmological perturbations. Since dynamical chaos is common in multi-field quar- 
tically coupled systems, I argue that the evolution of the fields after inflation must be examined 
to determine whether the amplitude of cosmological metric perturbations is underestimated in the 
standard inflationary calculations. 1 illustrate this effect with a simple hybrid inflation model. 
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1. Introduction. — A period of inflation has become the 
best candidate for the origin of structure in the universe. 
Recent years have seen considerable interest in the pro- 
cess of reheating, which follows inflation and transforms 
the inflaton field's large potential energy into particle ex- 
citations of various fields. A key development has been 
the understanding that parametric resonance can lead to 
extremely rapid particle production in regions of fc-space, 
a process known as preheating Q]. If the preheating dy- 
namics involves more than one scalar field, the resonance 
can extend to k = 0, so that field perturbations which 
are super-Hubble during preheating are amplified pl|§. 
Large scale metric curvature perturbations can be am- 
plified during preheating as well, which contradicts the 
usual assumption that they are constant on super-Hubble 
scales. The growth of super-Hubble metric perturbations 
is possible only when entropy modes are present — sin- 
gle field perturbations are purely adiabatic and hence the 
usual conservation law applies These large scales 

are relevant to structure formation and the Cosmic Mi- 
crowave Background (CMB), and hence it is crucial to 
understand this stage of the universe's evolution. 

In this Letter, I generalize previous studies of paramet- 
ric resonance in an attempt to understand under what 
conditions the growth of super-Hubble metric perturba- 
tions is possible. I point out a previously undiscussed 
general route for the amplification of super-Hubble scalar 
field and metric curvature perturbations in multi-field 
models, namely through dynamical chaos in the back- 
ground field evolution. (Note that "dynamical chaos" 
here refers to the evolution of a low-degree-of-freedom 
dynamical system of homogeneous fields, and should not 
be confused with "chaotic inflation" , which refers to the 
insensitivity of certain inflationary models to the inflaton 
initial conditions.) The idea is simple to state: if the ho- 
mogeneous background field dynamics is chaotic, then by 
definition small homogeneous perturbations of the back- 
grounds will grow exponentially on some characteristic 
time scale. This implies the growth of large scale metric 
perturbations, which are the modes relevant to structure 



formation and the CMB. Since dynamical chaos is com- 
mon in nonlinear systems with two or more degrees of 
freedom (e.g. in quartically coupled oscillators), I stress 
that the chaotic overproduction of super-Hubble modes 
may rule out certain inflationary models. 

It is important to point out the distinction between 
parametric resonance and dynamical chaos. While both 
can result in exponential production of super-Hubble 
modes, parametric resonance can be described analyt- 
ically while chaos cannot. Parametric resonance tech- 
niques can be applied in the special case that the back- 
grounds (or conformally scaled backgrounds) undergo pe- 
riodic motion. Essentially, linear stability analysis is 
performed about the periodic orbit, and the presence of 
instability implies exponential growth of perturbations. 
Dynamical chaos, on the other hand, involves instability 
to perturbations during extremely complex (and more 
general) evolution of the backgrounds. The lack of pe- 
riodic motion (and hence of parametric resonance) does 
not imply that perturbations cannot grow exponentially. 

There have been a number of studies of chaos in sys- 
tems of homogeneous fields in cosmology, although ap- 
parently none have made the connection with the growth 
of super-Hubble perturbations. Easther and Maeda 0] 
studied the chaotic dynamics of a two-field hybrid infla- 
tion system during reheating, although they did not in- 
clude metric perturbations. They found two effects: the 
enhancement of defect formation, and a significant vari- 
ation in the growth of the scale factor. However, they 
claim that only scales k ~ aH at the time of reheating 
will be affected. Cornish and Levin studied a single 
field model, and followed the evolution for several "cos- 
mic cycles" of bang and crunch. Chaos is possible in such 
a simple system if gravity is important, as we will see be- 
low. Latora and Bazeia || studied a class of two-field 
quartically coupled systems which are chaotic in some 
regions of parameter space. 

2. Perturbation dynamics. — A general dynamical sys- 
tem consists of a set of n (in general nonlinear) coupled 
first order ODEs 
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If we consider the state Xi to be a function of the initial 
state X® (and time), we can write 



5xi 



dxi{x®) 



5x» 



(2) 



In words, Eq. (||) says that once we have obtained a 
"background solution" Xi(x®) to the equations of mo- 
tion, we can obtain from it the evolution of a small per- 
turbation simply by taking derivatives. (Of course in 
specific cases, it may be impossible to obtain the back- 
ground solution in analytic form in order to implement 
this method.) 

The system (g) is said to exhibit dynamical chaos if 
the phase space is bounded and if a perturbation length 
d(t) = (Sx l Sxi) 1 ' 2 grows exponentially with time, i.e. 



d(t) ~ d(t )e 



hi 



(3) 



for sufficiently small initial displacement d(to) and suffi- 
ciently late t. Here h is known as the (largest) Lyapunov 
exponent || . The requirement of a bounded phase space 
excludes trivially unstable systems, such as the inverted 
harmonic oscillator. Necessary conditions for dynami- 
cal chaos in the system (0) are that it contain nonlinear 
terms and that n > 3 |9(], so that even very low dimen- 
sional systems can exhibit amplification of perturbations 
and the consequent extremely erratic background evolu- 
tion. 

I now consider the case of a spatially flat Friedmann- 
Robertson- Walker background universe with N mini- 
mally coupled real scalar fields 4>i described by the La- 
grangian density 



C 



(4) 



The equations of motion for the homogeneous back- 
ground scalar fields are 



3H<fn + V t = 0, 



(5) 



where Vi = dV/d(f>i. The background metric evolution 
is specified by the 0-0 Einstein equation, 



H 2 = 



8tt 
3m|j 



(6) 



Note that on the assumption that H does not change 
sign (which holds in realistic models where H > 0), we 
can substitute the square root of Eq. (|J) into Eq. (||). 
The evolution equations for the set of functions {</>j, 
can then be written in the form of the dynamical system 
(Q), with n = 2N. Thus with two or more nonlinearly 
coupled scalar fields, dynamical chaos is possible in the 



homogeneous background system. (Indeed we can see 
that if we allow H to change sign, as in bang-crunch 
scenarios, then we cannot eliminate H from Eq. (^). In 
this case n = 2N + 1, and a single scalar field is sufficient 
for chaos, as was found in M.) 

We can write the metric in the presence of scalar per- 
turbations in the general form pfl 



ds 2 = (1 + 2A)dt 2 - a 2 (t){2B ]i dx i dt 
-[(1 - 2ip)6ij + 2E\ lj ]dx l dx j }. 



(7) 



Here A, B, ip, and E are scalar functions, and subscript \i 
refers to the covariant derivative on the background con- 
stant time hypersurface. Two of the four scalar functions 
may be determined by a choice of gauge, and the linear 
perturbation evolution equations will reflect that choice. 
To obtain the long- wavelength limit of these equations, 
Sasaki and Tanaka pj (see also jl^]) noticed that it was 
sufficient to use Eq. (H) (in a slightly different form) . The 
issue of which gauge the perturbation equations are writ- 
ten in is tied to the choice of time variable used in the 
background equations ||llf . Explicitly, they found that to 
obtain the perturbed equations in a particular gauge, it 
is necessary to use a time variable which is not perturbed 
in that gauge. 

For example, in the synchronous gauge, A = B = 0, 
the usual cosmological time t itself remains unperturbed. 
Thus Eq. (||) applied to the background equations ^) 
gives 

5& + + 3H5'<j> + V tj 8<j) = 0. 

With the identification 



(8) 



SH = 4' 



(9) 



we obtain the usual perturbation equation of motion in 
the synchronous gauge, in the limit k — > 0. 

The importance of this result is that, if the homo- 
geneous background dynamics is chaotic, the quantity 
8<pi calculated above (which represents the homogeneous 
field perturbation in a particular gauge) will generically 
grow exponentially with time, as Eq. (|3j) indicates. Since 
spatial derivative terms are negligible on super-Hubble 
scales, this implies the exponential growth of super- 
Hubble finite wavelength modes. Finally, we can write 
the comoving curvature perturbation 1Z by H 



K = ip+ —6a, 
a 



(10) 



where a and 5a are the field and perturbation compo- 
nents along the direction of the background trajectory. 
Since the total length of the perturbation vector grows 
exponentially (unless the initial perturbation happens to 
be precisely adiabatic, corresponding to a perturbation 
in time), so must the component Sa, and hence so gener- 
ically must 1Z. 



2 



3. Hybrid inflation. — I can illustrate these ideas with 
a double scalar field model, which, as discussed above, is 
sufficiently complex for dynamical chaos to be possible. 
A popular class of two-field inflationary models is hybrid 
inflation In these models inflation can be termi- 

nated by a symmetry breaking transition in one of the 
fields, and the subsequent oscillations can be chaotic 0. 
1 will consider the potential 



\ U, \ ) = ^(Af 2 - A X 2 ) 2 + imV + i 9 W. (11) 



Inflation occurs at large (j>, where the effective mass of the 
X field, m 2 = g 2 (f) 2 + X\ 2 — M 2 , is large and the x field sits 
at the bottom of a potential valley at x = 0. Once the 
inflaton field drops below the critical value <f> c = M/g, 
the mass-squared m 2 becomes negative (the potential 
valley becomes a ridge), and the fields undergo a sym- 
metry breaking transition to one of the global minima 
at 4> = 0, X = where xo = M/vA. I consider the 

"vacuum-dominated" regime, where the potential during 
the inflationary stage, V{4>) = M 4 /(4A) + m 2 2 /2, is 
dominated by the false vacuum energy term M 4 /(4A). I 
also consider the case where the Hubble parameter at the 
critical point, 



H 2 



2tt M 4 

3A TTlpj 



(12) 



is much smaller than the oscillation frequencies about 
the global minima, which are = gM/ vA and rn x = 
V2M for small oscillations. This ensures that the fields 
will oscillate very many times after the critical point is 
reached before Hubble damping is significant. 

Preheating has been studied in hybrid models for vari- 
ous parameter regimes in the absence of metric perturba- 
tions |l5ll . The behaviour of large-scale metric perturba- 
tions was studied in |l6| , where it was found that growth 
is possible on large scales. Oscillations in hybrid inflation 
were found to be chaotic in although for a very dif- 
ferent parameter regime than I examine here. Also, the 
connection with the growth of large-scale perturbations 
was not made in 

Since my interest in this Letter is to establish a kine- 
matical connection between dynamical chaos in the back- 
ground fields and exponential growth of metric perturba- 
tions, I ignored the evolution of perturbations during the 
inflationary stage. It is important to note that for the 
parameters I consider, the large x mass during inflation 
implies damping of large scale perturbations during in- 
flation. Thus the question of whether the amplitude of 
metric perturbations produced in this model is consis- 
tent with the Cosmic Background Explorer normaliza- 
tion is not addressed here. A careful analysis, following 
the evolution of all important scales during inflation and 
preheating, and including the effects of backreaction, is 
required H. 



I considered the slice through parameter space speci- 
fied by M = 10~ 8 topi, m = 10 -16 m P i, A = 1(T 3 , and 
g 2 = 10 -2 -10 -4 . These parameters give an amplitude 
of cosmological density perturbations of the order 10 -5 
according to the standard inflationary calculation [ p^[ . 
I evolved the homogeneous background fields according 
to Eqs. (||) with initial conditions 4>(to) — 0.999</> c and 
x(*o) = O.OOlxo- I followed the evolution of the curvature 
perturbation on uniform-density hypersurfaces £fc (which 
coincides with 1Z on large scales) for k/a = lO~ 3 Ho using 
the longitudinal gauge equations 



k 2 
+ —t 



V 



4-7T 

— 2~ ( 

m Pl 



4$0, - 2V,&, (13) 



(14) 



Cfe = $fc - -r- ( §k + H$ k 
H 



(15) 



where $fc is the longitudinal gauge metric perturbation. 
I calculated the largest Lyapunov exponent h for the 
background evolution using Eq. (|§|). The results, shown 
in Fig. [|, indicate the presence of rich structure as g 2 
is varied. Regions of chaos with h ~ 0.05M are inter- 
spersed with regular stability bands where h ~ 0. The 
most prominent stability band is near the supersymmet- 
ric point <7 2 /A = 2. The correlation between the Lya- 
punov exponent and the growth rate of large scale metric 
perturbations is strong numerical evidence in support of 
my arguments. 

The growth at large scales is not simply due to the neg- 
ative mass-squared (or "tachyonic") instability near the 
potential ridge at x = [0. To demonstrate this, I plot 
in Fig. ^ the logarithmic growth rate of Sxk perturba- 
tions as a function of scale k. The solid line corresponds 
to the same parameters as in Fig. |l] (with g 2 = 10 -3 ), 
and shows a growth rate which is large at scales k ~ aM 
(due to the tachyonic instability) , and approaches a con- 
stant (the chaotic Lyapunov exponent) as k/{aH) — > 0. 
The dotted line shows the results for the same parameters 
except with the field trajectory constrained to <j) = 0. In 
this case, there is no large-scale growth (and no chaos) 
as expected for the effectively one-dimensional dynam- 
ics, while there is strong growth (due to the tachyonic 
instability) at small scales. Thus the two effects are dis- 
tinct, since both trajectories pass through the tachyonic 
instability region, while only one exhibits growth on large 
scales. 

4- Summary. — I have demonstrated a general link 
between instability in scalar field background evolution 
and the growth of super-Hubble metric perturbations. 
In particular, dynamical chaos in the fields can drive the 
growth — it is not necessary to have periodic motion and 
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FIG. 1. Largest Lyapunov exponent (solid line) for the ho- 
mogeneous fields, and logarithmic growth rate of (dotted 
line) for a scale k/a = W~ 3 Ho- The homogeneous fields 
are oscillating about one of the hybrid model's global min- 
ima. The parameters are M = 10~ 8 mpi, m = 10~ 16 mpi, and 
A = 1(T 3 . There is a clear correlation between the two curves. 
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FIG. 2. Logarithmic growth rate of 8\k as a function of 
scale for the parameters of Fig. |l| and g 2 = 10 -3 (solid line) 
and for the trajectory constrained to <f> = (dotted line). 



parametric resonance in the fields. Since chaos is com- 
mon in multi-field systems, it is important to examine 
the super-Hubble evolution during preheating carefully, 
and also to follow the evolution during the inflationary 
period, in order to determine whether the model conflicts 
with CMB measurements. 
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